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ABSTRACT 

Measuring the 3D distribution of mass on galaxy cluster scales is a crucial test of the ACDM 
model, providing constraints on the nature of dark matter. Recent work investigating mass dis- 
tributions of individual galaxy clusters (e.g. Abell 1689) using weak and strong gravitational 
lensing has revealed potential inconsistencies between the predictions of structure formation 
models relating halo mass to concentration and those relationships as measured in massive 
clusters. However, such analyses employ simple spherical halo models while a growing body 
of work indicates that triaxial 3D halo structure is both common and important in parame- 
ter estimates. We recently introduced a Markov Chain Monte Carlo (MCMC) method to fit 
fully triaxial models to weak lensing data that gives parameter and error estimates that fully 
incorporate the true shape uncertainty present in nature. In this paper we apply that method 
to weak lensing data obtained with the ESO/MPG Wide-Field Imager for galaxy clusters 
A1689, A1835, and A2204, under a range of Bayesian priors derived from theory and from 
independent X-ray and strong lensing observations. For Abell 1689, using a simple strong 
lensing prior we find marginalized mean parameter values Af2oo = (0.83 ± 0.16) x 10^^ 
H^^Mq and C = 12.2 ± 6.7, which are marginally consistent with the mass-concentration 
relation predicted in ACDM. With the same strong lensing prior we find for Abell 1835 
M200 = (0.67 ± 0.22) X 10^'" H-^Mq and C = 7.1^7"^^, and using weak lensing informa- 
tion alone find for Abell 2204 M200 = (0.50 ± 0.19) x lO^^ H^^Mq and C = 7.1 ± 6.2. 
The large error contours that accompany our triaxial parameter estimates more accurately 
represent the true extent of our limited knowledge of the structure of galaxy cluster lenses, 
and make clear the importance of combining many constraints from other theoretical, lensing 
(strong, flexion), or other observational (X-ray, SZ, dynamical) data to confidently measure 
cluster mass profiles. If we assume CDM to be correct and apply a mass-concentration prior 
derived from CDM sti-ucture formation simulations, we find {M200 — (0.99 ± 0.18) x 10^^ 
H-^Mq- C = 7.7 ± 2.1}, {Afaoo = (0.68 ± 0.19) x lO^^ /i-^M©; C = 4.4 ± 1.6}, and 
{M200 = (0.45 ± 0.13) X 10^5 /i-^M©; C = 5.0 ± 1.7} for A1689, A1835, and A2204 
respectively. 

Key words: gravitational lensing - cosmology: observations - dark matter - galaxiesxlusters: 
individual - methodsistatistical. 



1 INTRODUCTION 

Galaxy clusters are ideal laboratories in which to study dark mat- 
ter, being the most massive bound structures in the universe and 
dominated by their dark matter component (~ 90%). Constrain- 
ing the clustering properties of dark matter is crucial for refining 
structure formation models that predict bo th the shapes of dar k 
matter halos and t h eir mass funct i on (e .g. Navarro et al.J ( Il997h : 
lEvrard et alj i l2002h : ISahcall et all ( |2003|) : .Dahle ( .2006.) ). Several 
methods are used to measure galaxy cluster dark matter profile 
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shapes and halo masses on a range of scales, including X-ray and 
Sunyaev-Zeldovich (SZ) studies, dynamical analyses, and gravi- 
tational lensing. However, all of these methods typically require 
simplifying assumptions to be made regarding the shape and/or dy- 
namical state of the cluster in order to derive meaningful constraints 
from available data. Crucially, while most parametric methods typi- 
cally assume spherical symmetry of the halo, observed galaxy clus- 
ters often exhibit significant projected ellipticity and ha los in CDM 
structure formation simulations (e.g. 'Bett et al.l ( l2007t) (using the 
Millennium simulation); Shaw et al. (200q)) show significant tri- 
axiality in cluster-scale halos, with axis ratios between minor and 
major axes as small as 0.3. Understanding and accurately incorpo- 
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Figure 1. 34' x 34' R-band wide-field images of Abell 1689, 1835, and 2204, from the ESO/MPG WFI at La Silla. 



rating the impact of this on cluster mass and parameter estimates is 
crucial for accurate comparisons between measured cluster prop- 
erties and model predictions (e.g. the ACDM mass-concentration 
relation). 

Efforts to understand the impact of triaxiality in gravitational 
lensing and its potentia l role in exp l aining apparent discrepancies 
with CDM began when lOguri et al.l ^2005 ) applied a fully triaxial 
NFW model to the shear map of Abell 1689 to find that it is con- 
sisten t with 6% of cluster-scale halos. These continued as lGavazzil 
( 12005 ) showed that a triaxial NFW can reconcile parameter values 
derived from observations of the cluster MS2 137-23 to predictions 
from N-body simulations. 

More generally, ICorless & Kin3 j2007h demonstrated in the 
weak lensing regime that neglecting halo triaxiality in parame- 
terised fits of NFW models to dark matter halos with axis ratios 
significantly less than one can lead to over- and underestimates of 
up to 50% and a factor of 2 in halo mass and concentration, respec- 
tively. While extreme cases of triaxiality are rare, such halos can 
be much more efficient lenses than their more spherical counter- 
parts, especially when in confi gurations that hide most of the triax- 
ial shape along the line of sight. lOguri & Blandford j2008l) recently 
quantified this expectation, showing that the strongest lenses in the 
universe are expected to be a highly biased population, strongly 
favouring orientations along the line of sight, high levels of triaxi- 
ality, and apparent sphericity in projection; these biases correspond 
to a population of very triaxial prolate ("cigar"-shaped) halos ori- 
ented with their long axes approximately aligned along the line-of- 
sight. Thus, for strong lensing clusters - even those that are circular 
in projection on the sky such as Abell 1689 - spherical symmetry 
is an unjustified and error-inducing assumption. 

Further, even halos with less extreme axis ratios are inaccu- 
rately fit by spherical models. We expect the vast majority of galaxy 
cluster scale dark matter halos to be triaxial, and so it is important to 
includ e this expectation in model fitting. Moreover. lCorless & Kind 
( I2OO8 ) (herein CK08) showed that triaxiality leaves no quantifiable 
signature in lensing data (i.e. in the maximum-likelihood values 
obtained fitting models), and so it is generally not possible to deter- 
mine a priori the expected impact of triaxiality on mass and con- 
centration estimates for a given lensing halo. 

In the past, triaxial models have not generally been fit to weak 
lensing data because they cannot easily be well-constrained, largely 
due to the intrinsic limitations of lensing. Because lensing is deter- 
mined only by the projected mass density and shear of the underly- 
ing mass distribution, it is inherently impossible to fully constrain 
a 3D structure without imposing strong priors on the shape of the 



Table 1. Properties of the three lensing clusters 
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halo or supplementing lensing data with other data types more sen- 
sitive to line-of-sight halo structure, such as dynamical irrformation. 

Despite these difficulties, that triaxiality has been convinc- 
ingly demonstrated to be an important factor in model parameter 
estimation in lensing analyses make it important to directly include 
it in NFW fits to galaxy clusters. For studies of individual clusters 
doing so is crucial if claims regarding the validity of the ACDM 
paradigm based on NFW parameter estimates are to be meaning- 
fully evaluated. It is also of importance across populations to obtain 
more accurate distributions of galaxy cluster parameters, and espe- 
cially in measuring the galaxy cluster mass function and constrain- 
ing the scatter in mass-observable re lations. If the mass function is 
significantly sloped as expected (e.g. lEvrard et 'Zl( l2002h ). even the 
best-case scenario of a symmetric scatter of mass estimates due to 
neglected triaxiality would lead to an asymmetric shift in the calcu- 
lated mass function, because there are more low mass halos to shift 
up in mass than there are high mass halos to shift down (see e.g. 
ICorless & Kin3 (f2008b)). 

In CK08 we introduced a Bayesian Markov Chain Monte 
Carlo (MCMC) method for fitting triaxial NFW models to weak 
gravitational lensing data. We demonstrated that it accurately re- 
covers mean parameter values and errors across a population of tri- 
axial dark matter halos similar to those seen in CDM structure for- 
mation simulations. In this paper we make the first application of 
that triaxial fitting method to the weak lensing signals of three well- 
known galaxy clusters: Abell 1689, Abell 1835, and Abell 2204. 
We reanalyze preexisting lensing data obtained with 34x34 arcmin 
observations of the three cluster fields with the ESO/MPG Wide 
Field Ima ger, first presented and m o delled with spheric a l mas s 



models in 



Clowe & Schneide: 



li (l200ll) . 
^ ^20021) . 



IClowe & Schneideil ( l2002h . 

and lKing. Clowe & Schneideil 

The MCMC method allows for the simple and statistically 
robust imposition of prior probability function on the parameters 
of the triaxial model. Though some would rather avoid priors al- 
together, simple models include hidden priors stronger than those 
we employ in this work; for example, a "simple" spherical model 
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includes implicit (5-function priors on the triaxial axis ratios. Pri- 
ors are of such importance in this problem because models of 3D 
structures fit using only 2D lensing data are intrinsically undercon- 
strained. CK08 showed that carefully chosen priors are crucial to 
obtaining statistically accurate results across a population, but that 
slight offsets between the true distribution of parameters in nature 
and the distributions employed as priors lead only to a very small 
decline in the accurate performance of the triaxial fitting method. 
We therefore continue in this paper the investigation of the be- 
haviour and best use of various prior probability functions in fitting 
the weak lensing signals of galaxy clusters, in specific application 
to A1689, A1835, and A2204. 

An additional advantage of the MCMC method is that it al- 
lows for the straightforward and statistically-robust inclusion of 
additional constraints from other independent and complementary 
observations such as strong lensing, SZ, X-ray and spectroscopic 
studies. To demonstrate that capability, we investigate the inclusion 
of constraints from strong lensing and X-ray observations of the 
three clusters, imposing the external constraints as prior probabil- 
ity functions on the cluster model. 

Section 2 briefly describes the observational data, reviews 
lensing by triaxial dark matter halos and the MCMC triaxial fitting 
method, and discusses the various priors employed in this paper. 
Section 3 presents our results, and Section 4 concludes with some 
discussion. Throughout we assume a concordance cosmology with 
fim = 0.3, h = 0.7, and a cosmological constant f^A — 0.7. 



2 WEAK LENSING ANALYSIS 
2.1 Observations 

i?-band images of Abell 1689, Abell 1835, and Abell 2204, ob- 
served on May 29-30 2000 with the MPG/ESO Wide-field Im- 
ager are shown in Fig. [T] Table [T] lists the positions, redshifts, 
and background source densities used in the lensing analysis for 
each of the three clusters. Details of the observations, shape mea- 
su rement techniques, and bac kground galaxy selection are given 
in IClowe & Schneider ( 2001) (herein CSOl) for Abell 1689 and 
in Iciowe & Schneider (2002) (herein CS02) for Abell 1835 and 
Abell 2204. Since those publications, new photometric data were 
acquired from CFHT observations of Abell 1689 and Abell 1835 
dBardeau et ^ ( l2005h ). supplementing the original i?-band obser- 
vations with photometry in B- and /-bands for A1689 and V— 
and /- bands for A 1835. Abell 2204 was already observed in B 
and V with the WFI in CS02. This photometric information al- 
lows for colour-based cuts to be used in selecting the background 
galaxy population, in addition to the simple magnitude cuts used 
for A1689 and A1835 in CS02 and CSOl; the colour criteria used 
fo r the backgrou nd galaxies in the field of Abell 1689 are described 
in lClowd ( l2003h (herein C03). This allows for a much improved se- 
lection of background galaxies; most importantly, it allows the re- 
moval of faint cluster members, which are indistinguishable from 
faint background sources in single-band observations. These cuts 
were made very conservatively, at a level to exclude all galaxies 
with z < 0.5. For Abell 1689 all galaxies consistent within 2a of 
being at 2 < 0.5 were rejected; for Abell 2204 and Abell 1835 
only galaxies within la were discarded. This slightly less stringent 
standard was required in order to maintain an adequate background 
source density because the colour fields for those fields were shal- 
lower. The mean lensing redshift of the background galaxy popu- 
lations is calculated using photo-z information from the COSMOS 



field dMobasher etal] ilOOij) ) by making redshift and magnitude 
cuts on that catalogue identical to those applied to our cluster lens- 
ing catalogues. It is computed to be < Zs >— 0.8 for all three 
clusters, and this value is used throughout this paper. Additionally, 
the A1689 le nsing catalogue is u pdated to include the KSB PSF 
corrections of iBacon et al.l ( 12001 ), improving the accuracy of the 
shear measurements. These corrections were already in place in the 
CS02 analysis of the other two clusters. 

The weak lensing data are taken from an aperture 1.6' < r < 
18.0' in order to exclude the strong lensing regions at the centres 
of the clusters where k ~ 1. 

2.2 Weak Lensing by Triaxial Dark Matter Halos 

Weak lensing distorts the shapes and number densities of back- 
ground galaxies. The shape and orientation of a background galaxy 
can be described by a complex ellipticity e" , with modulus | = 
(1 — b/a)/{l + b/a), where b/a is the minor:major axis ratio, 
and a phase that is twice the position angle 0, e" — |e''|e^"^. The 
galaxy's shape is distorted by the weak lensing complex reduced 
shear, 17 = 7/(1 — ^), where 7 is the lensing shear and «: the con- 
vergence, such that the ellipticity of the lensed galaxy e becomes 



+ 7 



(1) 



l+3*e= 

in the limit of weak deflections and where * denotes complex con- 
jugation. Assuming a zero-mean ellipticity for the unlensed popu- 
lation, the expectation values for the lensed ellipticity on a small 
piece of sky is < e >= g « 7. This is the basis for weak lensing 
analysis in which the shapes of images are measured to estimate the 
shear profile generated by an astrono mical lens; a thoro ugh descrip- 
tion of weak lensing is given in Schne ider. King. & Er ben (2000). 

We calculate the 2D dispersion in the source ellipticities a^ — 
\J (jji + (T^2 for sach cluster catalogue from the measured ID dis- 
persion in the radial component of the background galaxy ellip- 
ticities measured after all colour and magnitude cuts, to which 
the gravitational shearing from the cluster contributes negligibly. It 
contains both the uncertainties due to the intrinsic noise in galaxy 
ellipticities and the measurement errors due to sky noise and PSF 
corrections. The values for each cluster are give in Table[T] 

A full parame terisation for a triaxial NFW halo is given by 
lling & Sutd ( I2OO2I) (herein JS02). They generalise the spherical 
NFW profile to obtain a density profile 

&cpc{z 



p{R) = 



(2) 



R/R,{l + R/RsY 

where 5c is the characteristic overdensity of the halo, pc the critical 
density of the Universe at the redshift z of the cluster, Rs a scale 
radius, R a triaxial radius 

y2 r72 

a'' 0-^ 

and a/c and b/c the minormajor and intermediate:major axis ra- 
tios, respectively. In a different choice from IS02 we define a triax- 
ial virial radius i?200 such that the mean density contained within 
an ellipsoid of semi-major axis -R200 is 200pc such that the con- 
centration is 
R200 



C 



Rs 



the characteristic overdensity is 
200 



Sr. = 



3 log(l + C)-^ 



(4) 



(5) 
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Figure 2. Prior probability distributions of axis ratios (Shaw), virial mass (Mass), mass-concentration (C), and triaxial line-of-siglit orientation angle (d), all 
derived from predictions of ACDM structure formation or lensing simulations. 



the same as for a spherical NFW profile, and the virial mass is 

M200 = ^^^abRloopc. (6) 

A more detailed de scription of this triaxial parameterisation and its 
benefits is given in lCorless & KinsI {2007). 

To allow comparison to some past lensing and numerical stud- 
ies employing spherical models, we define an effective spherical 
virial mass and concentration which we give in addition to the 
fully triaxial quantities when quoting marginalized and maximum- 
likelihood parameters. The effective spherical virial radius is de- 
fined as the radius r2oo at which the mean density within a sphere 
of that radius is 200 times the critical density, and the effective 
spherical virial mass the mass within that sphere: 

m2oo = (8007r/3)r2ooPc. (7) 

We further define the effective spherical concentration Csph of a 
triaxial halo as the ratio of the effective spherical virial radius to 
the geometric mean of the triaxial scale radii rs = Rs (abc)^''^: 

Csph = r2oo/rs- (8) 



Throughout this work we always directly fit the triaxial mass and 
concentration, then calculate the effective spherical values from 
those fitted triaxial models. For a spherical model, the triaxial and 
effective spherical parameterizations are identical: 771200 = M200 
and Csph = C. We consistently quote masses and concentrations 
derived using spherical models in our own work and that of other 
authors as the more general A/200 and C. 

The triaxial halo is oriented with respect to the observer by 
angles 9, the orientation angle between the major axis of the halo 
and the observer's line-of-sight, and 0; randomly oriented halos are 
distributed uniformly in </> and sin6'. 

The lensing properties of a triaxial halo, expressed as 7 
and K, are determined by performing integrals over the spher- 
ical NFW convergence a nd its derivatives. They are derived in 
lOguri. Lee. & Sutg ( |2003|) (h erein OLS), and we su mmarised and 
extended some of that work in ICorless & Ki"n3 ( l2007h : it is included 
here in Appendix lAl 

2.3 MCMC Fitting Method 

A more detailed description of the MCMC method for fitting tri- 
axial NFW halos is given in CK08. MCMC methods employ a 
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Table 2. The marginalized mean parameter values obtained fitting single triaxial NFWs to A1689, A1835, and A2204 under a range of prior probability 
distributions on the halo mass, concentration, shape and orientation derived from CDM structure formation simulations and external observational constraints, 
described in Sec. |2.4| la eiTors on the estimate of the mean parameter values are given; all values are for a concordance cosmology of Q^n = 0.3, Ha = 0.7, 
h = 0.7. 
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0.82 ±0.14 


1 90+0.30 
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"guided" random walk that returns a sample of points represen- 
tative of the posterior probability distribution; the probabihty of a 
certain region of parameter space containing the true model is di- 
rectly proportional to the density of points sampled in that region. 
From the distribution of sample points the full posterior probability 
distribution is obtained, which is easily and directly marginalized 
over to obtain fully marginalized mean most-probable parameter 
estimates for all parameters. Such methods are particularly valuable 
in under-constrained or highly-degenerate fitting problems such as 
this one of fitting 3D triaxial models to 2D lensing data. By con- 
trast, with such weak constraints a maximum likelihood approach is 
both impractical and of very limited scientific value; the posterior 
distribution will be very flat with significant degeneracies giving 
poorly constrained maximum likelihood values. By exploring the 
full posterior probability distribution, this MCMC method allows 
the derivation of parameter estimates (and their accompanying er- 
rors) that account for the true uncertainties when fitting parametric 
models to lensing data. 

To implement the method, we must define the posterior prob- 
ability function, defined in Bayesian statistics as 

p{e\Tv)p{Tv) 



P(0) 



(9) 



where p{6\it) is the likelihood £ of the data given the model pa- 
rameters (the standard likelihood), p{n) is the prior probability 



distribution for the model parameters (e.g. a distribution of axis 
ratios drawn from simulations), and p{6) is a normalising factor 
called the evidence, of great value in comparing models of differ- 
ent classes and parameter types, but computationally expensive to 
calculate and unnecessary for the accurate exploration of the poste- 
rior distribution. We define the log-likel ihood function in the stan- 
dard m ann er for weak lensing follow ing lSchneider, King. & Erbenl 
( I2OOOI) and lKing & Schneidej i200lh 



lnp,ieMe,;n)). 



(10) 



where the reduced shear g is calculated using the triaxial conver- 
gence and shear of Eauations lA12l n-y is the number of background 
galaxies used in the analysis, and 11 is a six-element vector of the 
parameters defining the model: triaxial virial mass A/200 i concen- 
tration C, minor axis ratio a, intermediate axis ratio b, and two 
orientation angles 9 and (j). We note that while the angles are de- 
fined over a range < 9 < n and < < 27r, because of the 
elliptical nature of the projected density contours, they give rise to 
unique lensing profiles only over the range of < 61 < 7r/2 and 
< < TT. The prior probability distribution, already mentioned 
in the introduction, will be discussed in detail in the next section. 

In our MCMC sampler we use the covariance matrix of an 
early run to sample in an optimised basis aligned with the degen- 
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Table 3. The approximate maximum-likelihood parameter values obtained fitting single triaxial NFWs to A1689, A1835, and A2204 under a range of prior 
probability distributions on the halo mass, concentration, shape and orientation derived from CDM simulations and external observational constraints, described 
in Sec. 12.41 Icr (68%) errors on the maximum-likelihood parameter values are given; all values are for a concordance cosmology of f2m = 0.3, = 0.7, 
h = 0.7. 
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eracies of the posterior. We tune the step sizes of the sampler to 
achieve an average acceptance rate of 1/3 in each basis direction, 
run three independent MCMC chains, started at randomly chosen 
positions in parameter space, for each lens, and sample the distri- 
bution space until the standard var(chain mean)/mean(chain var) 
indicator is less than .15, strongly indicating chain convergence 
jLevyis & Bridld ( |2002|)). We utilise the G etDist package from the 
standard CosmoMC 1 Levyis & Bridld [20031 distribution to calcu- 
late convergence statistics, confidence contours, and marginalized 
and maximum-likelihood parameter estimates. We employ 60 bins 
in the Gaussian smoothing of the contours; choosing such a large 
number leads to noisier contours, with more small islands of con- 
fidence and irregularly shaped confidence regions, but is the level 
at which the size of the contours is no longer sensitive to the bin 
number. If a lower number is employed, the contours look much 
smoother but grow significantly larger leading to overestimates of 
the confidence regions. 

Throughout the paper we quote both marginalized and 
maximum-likelihood parameter values. It is important to note that 
the maximum-likelihood values returned by MCMC methods such 
as this do not correspond to the true likelihood maximum - no min- 
imization to find the exact peak of the likelihood distribution is 
attempted because in our six-parameter, degenerate triaxial model 
such peaks are poorly constrained. Instead, the quoted maximum 



likelihood values are the parameter values at the point with max- 
imum likelihood in the MCMC sample. These will correspond to 
values close to those at the peak of the likelihood distribution, and 
are included to demonstrate the offset that may occur between the 
maximum likelihood values and the marginalized mean parameter 
values. This offset occurs because there are significant asymmet- 
ric degeneracies between the triaxial parameters and the errors on 
those parameters are highly non-Gaussian. The two types of "best" 
parameter values reflect slightly different information about the 
lensing cluster: the maximum-likelihood value gives the one point 
at which the very best fit is attained, while the marginalized values 
give the region of parameter space where most probability resides. 
The marginalized statistics are more robust than the maximum- 
likelihood values, especially for an underconstrained problem such 
as this in which the degeneracies are significant and the errors large; 
this is reflected in the much larger error bars on the maximum- 
likelihood values. We will therefore treat the marginalized param- 
eter values as the "best" parameter estimates for our clusters, but 
quote maximum-likelihood values as well to provide a better sense 
of the one best model and to allow for more direct comparisons 
with previous studies which typically quote maximum-likelihood 
or minimum ^-squared parameter estimates and errors. The errors 
quoted for the marginalized mean parameter estimates are the la 
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Figure 3. The marginalized mean (top row) and maximum-likeliliood (bottom row) A/200 and C values for A 1689, A 1835. and A2204 under a flat triaxial 
prior as well as priors on halo mass, concentration, shape, and orientation derived from ACDM structure formation simulations. Error bars gives Icr en'ors on 
the mean and maximum-likelihood values, and the dotted lines gives the mass-concentration relation predicted in ACDM at la. 



standard deviation; those for the likelihood estimates are the one- 
sided 68% la limits. 



2.4 Prior Probability Functions 

As noted earlier, because lensing is a fundamentally 2D phe- 
nomenon, models of 3D structures fit using lensing data alone ~ 
whether strong, weak, or flexion - are fundamentally undercon- 
strained. It may therefore be of interest to further constrain the tri- 
axial models with prior probability functions derived from theo- 
retical predictions or complementary observations. Several of these 
functions were discussed in CK08; some of those are included here, 
and others in this paper are improvements or additions to the set 
presented in that work. Importantly, CK08 showed that the triax- 
ial fitting method applied to a population of triaxial halos using a 
relatively accurate prior on the distribution of triaxial axis ratios 
returns statistically accurate and robust estimates of parameter val- 
ues across that population. However, given our limited knowledge 
of the universe (thus our desire to model dark matter halos in the 
first place) the choice of prior is not universal or obvious, and the 
choice must be carefully made with respect to the particular prob- 
lem at hand. For example, when testing the mass-concentration 
relation that dark ma t ter halos are pred icted to exhibit in CDM 
jNavarro etell l ll997l) ; iNeto et al.l ilOOH) ). adopting a prior based 
on that relation would be inappropriate. However, in another prob- 
lem, when modelling a large population of dark matter halos to 



calculate a mass function, a loose mass-concentration prior may be 
an appropriate way to incorporate existing knowledge of the clus- 
ter and group population into the model. Similarly, a strict prior 
on triaxial axis ratios may be well-justified when applied across a 
large population, but may be a poor choice when looking at a single 
strong-lensing cluster likely to be an outlier in the triaxial distribu- 
tion. 

The prior probability functions examined in this paper are de- 
scribed below, and plotted in Figure |2] All priors but the most gen- 
eral are fundamentally grounded in CDM, whether directly from 
CDM structure formation simulations or indirectly from work em- 
ploying NFW profiles derived from those structure formation sim- 
ulations. 

• Flat: The Flat prior is a very general prior on the shape and 
orientation of triaxial halos. It is defined (slightly more broadly than 
in CK08) to have probability p — 1 when 

0.3 s; b s; 1.0 and 

0.3 s; a/b s; 1.0, (11) 

and p = elsewhere. It excludes models with extremely small 
axis ratios because they are not expected in nature, and their inclu- 
sion slows convergence of the MCMC chains significantly. It also 
includes a prior on the orientation angles of the triaxial halo corre- 
sponding to random orientation with respect to the line-of-sight: 

p{e) = 0.5 sin 0; 
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(12) 



This prior is the most basic employed in this paper, and it is applied 
in every case - all other priors are implemented on top of it. 

• Shaw: The Shaw prior is a CDM prior on the triaxial axis ra- 
tios, employing the axis ratio d istributions f ound i n the CDM struc- 
ture formation simulations of IShaw et alj ( l2006h . These distribu- 
tions are fit using polynomials; these are plotted in the first panel 
of Figure |2] and given explicitly in Appendix IB] 

• Mass: The Mass prior is a CDM prior on halo mass, taken 
from the mass function fit to the ACDM structure formation sim- 
ulations of lEvrard et alJ |2002). The differential number density of 
dark matter halos as a function of mass and redshift is given by 



n{M,z) = 



Apr, 



M 



-Q,//(M)exp [-|lricr"'(A-f) + (13) 



where ae// is the effective logarithmic slope, a'^{M) is the vari- 
ance of the density field smoothed on scales enclosing mass M at 
the mean density pm{z), and A, B, and e are free parameters fitted 
for a ACDM universe in Evrard et al. This prior favours lower mass 
halos, as seen in the second panel of Figure|2] 

• Concentration: The Concentration (C) prior is a CDM 
prior relating halo concentratio n to mass, t aken from the mass- 
concentration relation derived bv lNetoetai] ( l2007l) from the simu- 
lated groups and clusters in the Millennium simulation: 



C : 



Mo, 



5.26 



1 + z \ 1014/i-iMe 



with a log-normal scatter 



p(logC) 



1 



27r 



exp 



log C- < log C> 



(14) 



(15) 



where < log C > is calculated via Eq. [14] for a given M200 and 
the dispersion a is taken to be 0.09 for masses less than 10^^ Mq 
and 0.06 for masses greater than that threshold, in rough agreement 
with the results of Neto et al. ( 2007). Examples of this prior at dif- 
ferent values of M200 are plotted in the third panel of Fig.|2] 

• Mass + Concentration: The Mass-Concentration (MC) prior 
combines the predictions of CDM for mass and concentration. It 
employs the Evrard mass function of the Mass prior, and a looser 
form of the mass-concentration relation of the C prior, plotted as the 
"Wide" prior in the third panel of Fig.[2] The dispersion of the log- 
normal probability distribution of concentrations at a given mass 
(T, as defined for the C prior in Eq. [15] is taken to be 0.18. That 
value, compared to those derived from simulations and used in the 
C prior above, is twice that for lower mass halos and three times 
that of higher mass halos. This increased flexibility makes this a 
more general (though still CDM-based) physical prior, that favours 
generally CDM-like halos but does not overpower the preferences 
of observed lensing data. 

• Sphere: The Spherical prior sets both axis ratios to unity, a — 
6 = 1. It is included for comparison to the spherical fits typical of 
most lensing analyses. 



3 RESULTS 

Under each of these priors we fit the lensing signal of each of the 
three clusters with a single triaxial NEW halo using the MCMC 
method, taking lensing data from an aperture 1.6' < r < 18.0'. 
We treat all background galaxies as though they are located on a 



sheet at redshift < Zs >— 0.8; the low redshifts of the lensing clus- 
ters justify this simplifying assumption, as only for higher redshift 
lenses that are in the heart of the redshift distribution is the distri- 
bution of source redshifts important (Seitz & Schneider 1997). 

Tables [2] and [3] give the mean marginalized and maximum- 
likelihood parameter values with la errors for all three clusters un- 
der each of the priors listed in Section [Z4I The corresponding ef- 
fective spherical parameters are also given, as well as the and re- 
duced of the maximum-likelihood triaxial fit to the lensing data. 
Figure [3] plots the mean marginalized and maximum-likelihood 
A/200 and C values, and the dotted lines show the predicted CDM 
mass-concentration relation with la errors. Figure [4] gives the 1- 
and 2a (68 and 95 percent) A/2ao — C confidence contours for 
each cluster and prior, again with the CDM mass-concentration re- 
lation overplotted in dashed lines. Figure [5] gives the marginalized 
posterior probability distributions of all six triaxial parameters for 
each cluster and prior. 

Priors are clearly very important, as the parameter mean val- 
ues, contours, and ID probability distributions shift significantly 
depending on the prior applied. The impact of each prior differs be- 
tween the three clusters, which is excellent evidence that the priors 
interact with the lensing data rather than overpower it. Examining 
the marginalized and maximum-likelihood parameter estimates, the 
ID marginalized distributions, and the M200 — C contours in com- 
bination gives insight into the way the various priors function. 

• Shaw: In every case the Shaw prior significantly changes the 
posterior probability distributions for the axis ratios, as expected 
since in its absence the constraints on the axis ratios are very weak, 
as the projected ellipticity and orientation angle constrain two axis 
ratios and two 3D orientation angles. However, it is important to 
note the case of Abell 2204, where even with only a Flat prior the 
lensing data proves very capable of constraining the minor axis ra- 
tio (this signifies a large ellipticity on the sky), and the Shaw prior 
almost completely overrides the data when it is applied. Thus, while 
the Shaw prior was shown in CK08 to be quite good for determin- 
ing the statistics of a large population of dark matter halos, it may 
suppress real information in the lensing data when interpreting in- 
dividual clusters. The large shift in the posterior probability distri- 
bution of the axis ratios forced by the Shaw prior leads to changes 
in the preferred orientation angles - see for example that in Abell 
2204 line-of-sight orientation angles close to S = are no longer 
favoured, and that the ID distributions for (j) flatten considerably 
for all three clusters. The impacts on the mass and concentration 
estimates are relatively small, and generally toward lower values. 

• Mass: The Mass prior prefers lower masses; in every case this 
leads to a movement of the A/200 — C contours to lower masses 
and higher concentrations, where the higher concentrations serve 
to boost the lensing signal to counter the reduction in mass. In all 
cases the 6 distributions move towards more lensing efficient line- 
of-sight orientations atO = 0, serving also to counter the decreases 
in mass induced by the prior. 

• C: The concentration prior functions very differently for the 
different clusters. In Abell 1689, where the contours under the Flat 
prior strongly favour concentrations higher than expected in CDM, 
and contain the CDM mass-concentration contours only barely at 
2a, the imposition of the C-prior completely dominates the lensing 
signal and shifts the contours dramatically to lower concentration 
values. By contrast, in the case of A2204 the C prior leads to only a 
small decrease in the marginalized concentration and no change at 
all in the maximum likelihood parameters, as the contours under the 
Flat prior are already consistent with the mass-concentration rela- 
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tion. Abell 1835 falls between the other two; the Flat prior favours 
somewhat (theoretically) high values of C, and the imposition of 
the C prior leads to a small shift of the AI200 — C contours to 
slightly lower concentration values. In every case, the marginal- 
ized mass moves to counter the prior-induced shift in C, moving to 
higher (lower) masses as C decreases (increases). 

• MC: As expected, this combination of the mass and (looser) 
concentration priors combines the effects of its components: in ev- 
ery case the marginalized mean parameter values lie in-between 
those of the individual Mass and C priors, as seen in Fig. [3] The 
MC prior generally moves the parameter estimates of the Flat 
prior towards the CDM mass-concentration relation, but weakly: 
the M200 — C confidence contours under the MC prior are almost 
entirely contained within those of the Flat prior in all cases. Still, 
however, the shift in mass and concentration parameters is quite 
dramatic for the outlier Abell 1689. Also note that the ID distri- 
butions for the axis ratios a and b and orientation angles and (f} 
typically fall between those of the two components, but usually fol- 
low those of the Mass prior more closely. 

• Sphere: The spherical prior generally prefers lower masses 
and concentrations to the Flat triaxial prior; this is clearly seen in 
the A'hoo — C contours, as the Flat contours are enlarged from 
their spherical counterparts in all directions but preferentially to- 
wards higher masses and concentration. This is because there are 
more inefficient lensing orientations for a triaxial halo than there 
are efficient ones. Thus, as argued already in CK08, halo triaxiality 
cannot explain a whole population of apparently over-concentrated 
halos, though it can quite easily account for some individual cases. 
Under this prior the estimated error contours on the parameters are 
much smaller than those under the Flat or Mass priors, but not as 
constraining as the concentration priors C and MC. 

The varied results for each of the clusters under the differ- 
ent priors indicate that the priors are not generally dominating the 
lensing signal. The stringent Shaw and C priors are seen to some- 
times work strongly against the data, but this is less of a concern 
because the impacts are so obvious - for example, the Shaw prior 
may be valid for Abell 1689, for example, where the marginalized 
ID distributions for a and h are entirely enclosed by their counter- 
parts under the Flat prior, but it is clearly not for Abell 1835 and 
especially Abell 2204 where they sit mostly outside the more gen- 
eral distributions. Thus the Shaw prior is often not a good choice 



for studies of individual clusters; however, it may still be the best 
choice for statistical studies, as it was shown to perform very well 
across populations in CK08. Similarly, the C prior applied to Abell 
1689 shifts the contours dramatically outside of those of the Flat 
prior, indicating it is a bad choice in this problem. Though the Mass 
prior does not behave so badly, it is an unbalanced expression of the 
CDM model, favouring lower masses while pushing up concentra- 
tion. Generally, the MC prior is the best behaved, weakly reflecting 
true prior expectations for any galaxy cluster, exerting a balanced 
influence on both mass and concentration while never entirely over- 
riding the data. 

In the case of Abell 1689 we find evidence supporting a con- 
centration significantly higher than that predicted by CDM sim- 
ulations. Every prior but the C and MC priors, which explicitly 
forces the solution towards the predicted mass-concentration rela- 
tion, gives marginalized and maximum-likelihood values of con- 
centration well above the Icr CDM mass-concentration corridor. In 
the i\f200 — C contours the weak lensing data of A 1689 are seen to 
be at best marginally compatible with the CDM mass concentration 
relation at 2a for all but the C and MC priors. 

The reduced values are slightly greater than unity for all 
three clusters under all priors, indicating the NFW density profile is 
a good fit to the lensing data; the values are expected to be slightly 
larger than unity, as they are, because our models do not fit any 
substructure. 



3.1 Cluster-Specific Priors 

In addition to these general priors, if independent constraints on 
mass and concentration exist for a cluster from X-ray, strong lens- 
ing, or dynamical data, it may be desirable to combine them with 
the weak lensing analysis to derive joint constraints. In this sec- 
tion we examine what external data is available for each of these 
clusters and how it may be combined to tighten constraints on the 
triaxial NFW. Care is necessary in doing so, for in some cases the 
very assumptions of sphericity we are trying to avoid may be im- 
plicit in the parameter values derived from other methods, and so 
combination with these results may be detrimental to the accuracy 
of our final mass and concentration constraints. 
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Figure 5. The ID marginalized probability distributions for all six parameters of the triaxial NFW fit to A1689, A1835, and A2204 under a Flat triaxial prior 
as well as priors on halo mass, concentration, shape, and orientation derived from ACDM structure formation simulations. AU priors are described in Sec. |2.4| 
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Figure 6. The left-hand panel plots the 2D marginalized 68 and 95 % confidence contours in the C — M200 plane for A1689 under the Flat, 8, 6e, and MC 
priors, described in the text. The dotted lines show the region of parameter space consistent at Icr with the mass-concentration relation predicted in ACDM. 
The right-hand panels give the con'esponding ID marginaUzed probabiUty distributions for all six triaxial parameter under each prior. 
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Figure 7. The ID marginalized distribution for the triaxial line-of-sight ori- 
entation angle d For Abell 1689 under the general CDM MC prior alone 
(plotted in solid black) compared with that under the MC prior in com- 
bination with the observational Einstein radius 9^ prior (plotted in dashed 
red). The addition of the Einstein radius constraint to the general CDM 
p rior favors more line - of-sig ht oriented models, supporting the predictions 
of lOguri & Blandford lioog) . 



3.1.1 Abell 1689 

Abell 1689 presents a wealth of strong lensing features, and these 
have b een modelled ex t ensively (e.g. [M iralda-Escud e & Babull 
( Il995l) . iBroadhurst et"ZI l l2005ahlLimousin et al.l (120071) ). Though 
many studies publish NIW mass and concentration parameters, 
spherical symmetry is always assumed in their derivation, and so 
direct application of these constraints as a prior on our triaxial 
model is not sensible. However, a model independent quantity that 
is very reliably constrained by the strong lensing work is the Ein- 



stein radius 9e (the location of the lensing critical curves where 
giant arcs appear); for Abell 1689 is consistently measured to 
be 45" for a spectroscopically constrained mean source redshift 
Zs = 1. There are two possible ways to apply this constraint: 



• lOguri & Blandford ^20081) demonstrate that galaxy clusters 
with relatively large Einstein radii such as Abell 1689 are strongly 
biased towards line-of-sight orientations where 6 is close to zero. 
They express this bias as the fraction of halos with Einstein radii 
above a cutoff value for which | cos 9\ > 0.9: for the Einstein ra- 
dius of Abell 1689 this fraction is approximately 80%. We give 
this probability distribution analytic form with a Gaussian peaked 
at 6 = with a dispersion that gives the correct fraction of systems 
with I cos > 0.9: 



ag \/2tv 



exp 



I cos 01 - 1)^ 



(16) 



where ag = 0.115. This prior, effectively a CDM prior on the ori- 
entation of strong lensing halos, which we will call the prior, is 
plotted in the fourth panel of Figure[2] compared with the distribu- 
tion for unbiased random halo orientations. 

• Einstein Radius: Because the mean source redshift is well- 
determined for the Einstein radius of Abell 1689, it is also possible 
to more directly constrain the weak lensing models, demanding the 
models give an Einstein radius for Zj, = 1 of 45", allowing for 
a 10% error in the measurement of 8e due to uncertainties in the 
redshift of the lensed source or the impact of baryons at the very 
centre of the cluster (ae„ — O.IOe): 



p(n) 



1 



cg^ 



exp 



2ai 
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The 2D AI200 — C confidence contours and ID marginalized 
parameter distributions under both of these priors are shown in Fig- 
ure [6] as well as under both priors in combination; the resulting 
marginalized and maximum-likelihood parameter estimates are in 
Tables[2]and[3] 
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The application of the 6 prior significantly shifts the contours 
to lower values of mass and concentration; this is expected because 
it explicitly favors more efficient lensing orientations. It also no- 
ticeably increases the probability of more extreme values of the in- 
termediate axis ratio 6. In the line-of-sight orientations favored by 
the 6 prior, circular symmetry in projection is achieved as a and b 
become close to equal in value, so for the relatively circular A1689 
b will follow a more closely under this pri or. Lower values of the 
axis ratios are favored b e cause , as shown in ICorless&KinibOOl) 
and lOguri & Blandforj ( 1200 Sh . small axis ratios give the most ef- 
ficient lensing for halos in line-of-sight orientations. As expected, 
the distribution of 6 values shifts towards zero, and this tighter con- 
straint on one orientation leads to a corresponding tightening of the 
distribution of its counterpart (f). 

The Einstein radius 9e prior shifts the A/200 — C contours dif- 
ferently, towards slightly higher concentrations and lower masses. 
The move away from higher masses suggests that increasing con- 
centration is a more efficient way of boosting the Einstein radius 
than increasing the mass. This preference for higher concentrations 
as opposed to higher masses also is consistent with the fact that it is 
the center of the cluster that determines the strong lensing behavior 
as opposed to the overall profile. This is important because it is at 
the cluster center that the impact of baryons, whether in the cen- 
tral cluster galaxies or the intracluster medium, are expected to be 
most important. Thus, this more stringent Einstein radius prior may 
not be well-suited for application to the fitting of a dark-matter only 
profile, as it may enhance the impac t of bar yons at the dense cluster 
center . However, iMeneghetti et al.l ( l2003h and IWambsganss et all 
( l2008 f) have both shown the impact of the central cluster galaxy on 
the strong lensing cross section to be significantly scale dependent: 
it may increase the cross section by up to 100% for small separation 
systems, but has a negligible impact for large separation systems. 
Thus, for a system such as A 1689 with a very large Einstein ra- 
dius, the contribution of baryons to the size of the Einstein radius 
is predicted to be small. 

We also combine the two strong-lensing priors, applying the 
orientation bias 9 and the Einstein radius 9e prior simultaneously. 
As might be expected, the resulting M200 — C contours fall be- 
tween those of the two component priors, picking out the lowest 
mass and concentration halos that generate an Einstein radius as 
large as that observed for the cluster. This is the most physically 
representative of the strong lensing priors: it combines both our 
expectations for the orientation of the triaxial halo given its role 
as one of the strongest lensing clusters in the universe with the 
demand that the fitted models adequately reproduce the observed 
strong lensing signal. Under it, A1689 is not consistent with the 
CDM mass-concentration at over 2a. However, as with the 9e prior 
applied alone, the potential impact of baryons at the cluster center 
is not accounted for. 

As a further test, we combine the observational Einstein ra- 
dius prior 6e with the general CDM MC prior. The results com- 
plete a consistent picture: adding an additional CDM constraint 
on the mass and concentration to the observational 6e prior shifts 
the distribution of the orientation angles to strongly favor orienta- 
tions near the line-of-sight {9 small), thus recovering the alternate 
CDM prediction, as implemented in the 9 prior, of a strong orien- 
tation bias toward low 9 values. The observational Einstein radius 
prior increases the preference for low 9 values compared to that 
under the MC prior alone, as shown in Figure|7j this provides ob- 
servational support for the link between large Einstein radii and 
line-of-sight triaxial orien tations in a CDM universe predicted by 
lOguri & Blandforj j2008h . 



Under the 9 prior alone, the 2a Ahoo — C confidence contours 
enclose the predictions of CDM for the mass-concentration rela- 
tion. With the addition of the 9e prior, however, the 2(7 i\/2oo — C 
confidence contours pass far above the predicted CDM mass- 
concentration corridor; the combined 9e-\-AIC prior also excludes 
the predicted mass-concentration relation. Note that under the 9 
prior both the triaxial mass and concentration estimates, and un- 
der the combined prior the triaxial mass only, are lower than their 
counterparts under the Spherical prior; this indicates that use of the 
spherical model may indeed lead to false conclusions regarding the 
place of A 1689 within the CDM paradigm. Due to remaining un- 
certainty regarding the potential importance of baryons at the clus- 
ter center, especially if A1689 has undergone a recent merger as 
the dynamical and X-ray data suggest it may have, we hesitate to 
quote results employing the Einstein radius prior as the definitive 
triaxial model for the cluster. We therefore take as the best current 
constraint on the triaxial NEW model of A 1689 the values derived 
under the less-restrictive strong lensing orientation 9 prior alone. 
These are M200 = (1.18 ± 0.23) x 10^^ M© and C = 12.2 ± 6.7 
(Mes = (l-16lo.3o) X lO^'' Mo and C,ph = 12.1 ± 6.8), which 
are marginally consistent at 2a with the prediction of CDM. How- 
ever, the size and shape of the A/200 — C confidence contours that 
stretch towards very high values of concentration while only graz- 
ing the predicted mass-concentration corridor, combined with the 
complete exclusion of that corridor when the triaxial model is fit 
under the combined 9 + 9e and 9e + MC priors, indicates that tri- 
axiality alone does not fully explain A1689's unusual characteris- 
tics: comparison of it with the other two clusters in this study shows 
that A 1689 is still unusual for the high masses and concentrations 
its lensing data favor. 

Abell 1689 has also been observed in the X -r ay us - 
ing both Chandra and XMM-Newton iXus & Wul ( l2002h ; 
lAndersson & Madeiskil ( l2004l) ; lLemze et al.l ( l2008h ). However. X- 
ray mass models require the assumption of both spherical sym- 
metry and hydrostatic equilibrium. Generally, because the gravi- 
tational potential traced out by the X-ray gas is always rounder 
than the underlying density distribution. X-ray models assum- 
ing spherical symmetry are less likely to be strongly affected 
by triaxi ality and elongation of the cluster along the line-of- 
sight (see lKassiola&Kovneil ([1993); Russell et al. (2008)). How- 
ever, the strong orientation bias of lOguri & Blandforj ( l2008h in- 
dicates Abell 1689 is likely to be very elongated along the line- 
of-sight, making sphericity a particularly poor assumption in this 
case, even though the X-ray isophotes are quite round. Further, 
lAnderssoir& Madejs^ j2004h find evidence of either large bulk 
motions of the intracluster gas or of an infalling subcluster, either 
of which makes it very unlikely the cluster is relaxed and in a state 
of hydrostatic equilibrium. Because there are no X-ray halo mod- 
els that account for this observed substructure, elongation, and dis- 
equilibri um, we do not at t his time include an X-ray constraint. 
Though iLemze et al. I j2008l) have previously attempted to com- 
bine lensing and X-ray observations of Abell 1689, their analysis 
discards the X-ray temperature profile because it is irregular, fur- 
ther suggesting an unrelaxed state of the cluster, and neglects any 
non-spherical structure of the cluster. Thus, at this time it seems a 
consistent and meaningful combination of X-ray and lensing con- 
straints is not yet possible. 

In addition, iLokas et"al] ( |2006|) compiled a spectroscopic 
study of the cluster from existing observations, but it did not in- 
clude enough position-velocity pairs to constrain more than the 
simplest models. It did however further confirm the complex state 
of Abell 1689, showing evidence for at least one associated struc- 
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ture in addition t o the m ain halo (it should be noted, however, 
that iLemze et alj ( l2008bh very recently carried out a larger spec- 
troscopic study of the cluster, in which they also identify an in- 
falling subclump but argue that otherwise the cluster appears uni- 
form with little substructure). While one subclump could perhaps 
be accounted for by the inclusion of a second clump in the weak 
lensing mass model in addition to the primary triaxial halo, we 
found the error contours grow so large when a second spherical halo 
is added as to be physically insignificant. [Limousin et al.l i2007t) 
present strong evidence of a second clump in their strong lens- 
ing analysis of the cluster and tightly constrain its position (in the 
same region where the X-ray data also suggests a major substruc- 
ture); we attempt a two-component fit with a spherical clump at the 
location of their Clump 2 in addition to the primary triaxial halo 
with a 1' positional uncertainty; we find that even with the position 
of the subclump so tightly constrained the weak lensing data has 
great difficulty constraining the model parameters. After running 
the MCMC chains for over four times as long as for the single tri- 
axial halo, convergence is achieved according to the standard statis- 
tics, though the likelihood distributions remain very un-smooth, in- 
dicating poor constraints. With that caveat, the marginalized mean 
mass of the main halo is reduced to 0.92 x 10^^ Mq while the 
secondary clump acquires mass 0.29 x 10^^ Mq, with respective 
concentrations of C = 21.5 ± 19.4 and C = 12.8 ± 12.6. The er- 
rors on the main halo concentration increase by a factor of 2 from 
the Flat triaxial case, and the 68% confidence contours for the sec- 
ondary clump include M — 0. Thus, the weak lensing data alone 
cannot meaningfully constrain a subclump so close to the cluster 
centre, even when its position is externally constrained. 

The inclusion of a second clump is not crucial for the weak 
lensing model, as the secondary peak i s located about 1.4' from 
the cluster centre jLimousin et al.l OOOTt) ). inside the central region 
excluded from the weak lensing analysis. Though weak lensing is 
a non-localized phenomenon and thus would not be entirely unaf- 
fected by such a clump, the effects would be small compared to 
those in strong lensing or X-ray analyses that focus on the very 
centre of the cluster. Additionally, if the second clump is indica- 
tive of a merger along the line-of-sight, which the dynamical and 
X-ray data suggest it may be, we ech o the argument already put 
forward in Oguri & Blandford l l2008h that a very prolate triaxial 
halo is a good zeroth order approximation of such a system. In- 
deed, the reality must be somewhere between the extremes of two 
separate clumps and a fully merged cluster, so a model treating the 
cluster as single triaxial merged halo is as good or better than any 
other model currently available to weak lensing analysis, until the 
dynamical structure of Abell 1689 is better constrained. 

In addition to having at least one associated structure, Abell 
1689 is likely part of a large-scale matter filament of the cosmic 
web oriented along the line-of-sight. This would add significant 
substructure to the system that cannot be accounted for entirely by 
a single triaxial halo or the addition of a single associated structure, 
and may contribute to the apparently high concentration of the sys- 
tem. We discuss large-scale line-of-sight structure further in |4.3l 



3.1.2 Abell 1835 

Abell 1835 has four identified strong lensing features, but these are 
not spectroscopica lly confirmed (Schmidt, Allen, & Fabian (2001); 
ISmithetal.l ( l2005h ). These suggest an Einstein radius of 31.3" for 
a source redshift of anywhere from z — 0.6 to z = 3. With such 
poor constraints on the source redshift, application of a direct Ein- 
stein radius prior is not possible. However, because the orientation 



bias is relatively robust to source redshift - though it does depend 
on source redshift, the difference in the percentage of halos with 
I cos9\ > 0.9 for a given Einstein radius is only ~ 5%, smaller 
than the uncertainties in the bias from uncertainties in the cosmo- 
logical parameters - we apply the 6 prior as a valid constraint de- 
rived from the strong lensing features of the cluster For the Einstein 
radius of Abell 1835, approximately 60% of the halos should have 
I 008 6*1 > 0.9: we implement this as for Abell 1689 via Eall6l 
taking as = 0.22 to reflect the weaker orientation bias for Abell 
I835's smaller Einstein radius. 

The 2D A4200 — C confidence contours and ID marginalized 
parameter distributions under the 6 prior are shown in Figure [8j 
the resulting marginalized and maximum-likelihood parameter esti- 
mates are given in Tables|2]and[3] Again, as for Abell 1689, the con- 
tours shift to favour lower masses and concentrations, and smaller 
values of b are favoured while (j) is more tightly constrained. 

The core of Abell 1835 has been observed in the X- 
ray with Chandra dSchmidt. Allen & Fabiad ( 1200 ih ) and XMM 
jMaierowicz. Neumann & Reiprichl ( l2002h ). The orientation bias is 
weaker for this halo because of its smaller Einstein radius, and 
our weak-lensing-only triaxial fits indicate no particular prefer- 
ence toward extreme axis ratios, making the assumption of spher- 
ical symmetry of the gravitati onal potentia l in th e X-ray analy- 
sis less problematic. Moreover. ISmith et al.l ( |2005|) argue that the 
gas properties of the cluster indicate it is relaxed, and there is 
no significant substructure apparent in lensing or X-ray analy- 
ses. Therefore, the X-ray mass and concentration can reasonably 
be employed as valuable constraints on the triaxial halo models. 
To do so, we convert the best-fitting NEW parameters found in 
ISchmidt. Allen. & Fabiai] (12001) for an Einstein-De Sitter cosmol- 
ogy to our ACDM cosmology, such that their best-fit parameters 
{raoo = l-'28tof2 ^"'Mpc; C = 4.0l^:^l} become {raoo = 
l-39to:34 /i"^Mpc; C = S.dtoil}- From these best-fitting spher- 
ical parameters we calculate the X-ray determined virial mass of 
the cluster to be M200 = l-lSl^g x 10^^ Mq. (Though there was 
some discrepancy in the temperature profiles derived from Chandra 
and XMM due to a p otentially unaccou nted for solar flare in the 
Chandra observation dMarkevitchl (l2002h). the virial mass derived 
from the XMM analysis bv i Maierowicz. Neumann & Reipricl] 
(2002) is fully consistent with the Chandra mass at la, so we em- 
ploy only the NEW fit from Chandra rather than the unparame- 
terized constraint from the XMM analysis). We apply these con- 
straints on the cluster mass and concentration (not on R200, since 
the triaxial and spherical definitions of radius differ) to the triaxial 
fit to the weak lensing data. We assign the error distributions Gaus- 
sian forms, with errors twice those given for the parameters in order 
to allow for small residual systematic errors due to the assumptions 
of spherical symmetry and hydrostatic equilibrium. 

Figure [8] shows the A/200 — C confidence contours and ID 
marginalized parameter distributions under the X-ray prior; the re- 
sulting marginalized and maximum-likelihood parameter estimates 
are given in Tables |2] and [3] Even with doubled error bars the X- 
ray signficantly tightens the constraints on both the mass and the 
concentration. It very slightly shifts both a and b to higher values, 
suggesting the spherical assumption of the X-ray model is having 
some impact on the combined parameter estimates, but that it is 
small. Similarly, a shift towards smaller values of implies a slight 
preference for line-of-sight orientations in which the halo is more 
likely to look circular on the sky. 

The X-ray and weak lensing data agree very well. In par- 
ticular. X-ray analyses, with their focus on the central regions of 
the cluster, typically do very well at accurately constraining the 
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Figure 8. The left-hand panel plots the 2D marginalized 68 and 95 % confidence contours in the C — M200 plane for A1835 under the Flat, 0, and X-ray 
priors, described in the text. The dotted lines show the region of parameter space consistent at Icr with the mass-concentration relation predicted in ACDM. 
The right-hand panels give the corresponding ID marginalized probability distributions for all six triaxial parameter under each prior. 



concentration but have more p r oblem s accurately constraining the 
mass (see e.g. [Limousin et alj ( l2007l )): thus the fact that it is the 
triaxial concentration that is most affected by the application of the 
X-ray prior, while the triaxial mass stays much the same, further in- 
dicates that the X-ray and tensing data are in good agreement, and 
that the application of the X-ray prior is reasonable and improves 
the final parameter estimates. 

Figure[8]also plots the Af200 — C confidence contours and ID 
marginalized parameter distributions under a combined 6 and X-ray 
prior. The mass and concentration are both very well constrained, 
as the preference for line-of-sight orientations slightly lowers the 
mass and concentration estimates from those under the X-ray prior 
alone. 



3.1.3 Abell2204 

Abell 2204 presents no known strong lensing features, and so 
neither the 9 01 6e priors may be applied. It has been ob- 
served in the X-ray w ith Chandra (Reiprich & Boehringer (2002); 
ISanders, Fabian, & T aylor (2005)), XMM l lZhang et al J 120081) , and 
Suzaku iReiprich eTal. C2008.) ) where it is seen to exhibit complex 
morphology potentially indicative of a recent merger. The com- 
plex shape of Abell 2204 is further evidenced by our weak-lensing 
only fits of the triaxial model, which under almost all priors favour 
highly triaxial halos with extreme axis ratios. Abell 2204 is thus 
unrelaxed and of complex shape, making assumptions of sphericity 
and equilibrium highly suspect; we therefore forego combination 
of the weak lensing constraints with X-ray observations. 



4 DISCUSSION 

We first compare our fits under the Spherical prior to those of pre- 
vious work on the three clusters - that used weak lensing, strong 
lensing, and X-ray data - to ascertain the agreement of our lensing 
catalogue and analysis with other methods and observations. Once 



the accuracy of our method and catalogues is so established, we 
then move on to discuss the impact of triaxiality on the measure of 
the mass and concentration of A1689, A1835, and A2204. 

4.1 Previous worls 

Table |4] collects the NFW M200 and C values published in 
past work on the three clusters, all converted to a co ncor- 
dance cosmology with h = 0.7; [Limousin et al. I ( I2OO7I) and 
lUmetsu & Broadhursll Jioos !) have previously summarised the ex- 
isting work on Abell 1689 in a similar way. Figure [9] plots most 
of those results on top of the confidence contours for the fits un- 
der a Spherical prior in this work. Though the Spherical prior fits 
are NOT our best estimates of the cluster parameters, they are used 
here in order to more directly compare our work to existing work in 
order to establish the consistency of the method. Once this compar- 
ison establishes our results as consistent and reliable, we will move 
on to discuss the fully triaxial results. 

4.1.1 Abell 1689 

The concentratio n of Abell 1689 has b een a source of signifi- 
cant controversy; [Limousin et al. I j2007l) offers an excellent and 
thorough discussion of the systematic issues that may affect con- 
centration estimates in the case of this remarkable cluster. Cru- 
cially, contamination of lensing catalogues with very faint clus- 
ter members is expected to lead to underestimates of the concen- 
tration when magnitude cuts alone are used to distinguish back- 
ground from foreground and cluster galaxies. Such was the method 
used in the early weak lensing s tudies of Abell 1689 in CSOl and 
iKing, Clowe & Schneider! ( l2002h . A version of that same catalogue 
is employed here, similar to that used in C03, improved with two 
bands of additional colour data and small improvements to the 
KSB PSF corrections. The mean lensing redshift of the background 
galaxy population is changed in this analysis from < >= 1 to 
< Zs >— 0.8, to account for the more stringent colour cuts in 
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Figure 9. The black contours give the 68 and 95 percent confidence contours for fits under a Spherical prior to the three clusters from this work, compared 
with the results of other authors fitting spherical models. Circles plot weak lensing-only results, diamonds weak lensing IR results, squares combined weak 
& strong lensing analyses, stars strong lensing-only measurements, and crosses X-ray results. The results and authors are compiled in Table|4] The Spherical 
contours do not represent the best estimates from this work, but are given for more direct comparison with past results to establish the consistency of the 
MCMC method with previous work. 



the galaxy population and make use of the improved photo-z infor- 
mation now available from the COSMOS field. Applying the same 
fitting technique as in CSOl and C03 to the updated catalogue, we 
find maximum-likelihood parameter values ofM2oo = 1.47x10^^ 
Mq and C = 11.1. The higher mass value compared to the C03 fit 
can be attributed primarily to the change in the mean background 
source redshift assumed in the analysis: overestimating the mean 
background redshift as was done in C03 shifts the A/200 — C con- 
tours to lower masses. In our MCMC maximum-likelihood fit under 
a spherical prior, we obtain marginalized mean parameter values 
M200 = 1.27±0.24x 10^^ Mq andC = 13.3±7.6, and approx- 
imate maximum-likelihood parameters Af2oo = l-SSlg gl ^^^^ 
Mq and C — 10. 6^4 q which are consistent with this most recent 
result at la. The concentration values have increased with every 
new iteration of the catalogue, as expected as the level of contami- 
nation by faint cluster members decreases. 

iKing et al.l ( l2002bl) performed a weak lensing analysis em- 
ploying infrared data. While the infrared provides an order of mag- 
nitude fewer background galaxies, it allows for a more secure de- 
termination of the cluster red sequence and removal of foreground 
galaxies. The mass and concentration yielded in that study are 
both lower than those we find here; however, that is expected be- 
cause that analysis did not include the ~ 1.1 shear correction of 
lBaconetal.l ( l200ll) . 

Our marginalized s pherical mass agrees w ell with those of the 
weak lensing results of ISardeau et al.l KOOSi) and [Limousin et ahl 
( l2007h . both using CFHT observations of the cluster, and is gener- 
ally lower than those using the Subaru observation. However, our 
mean concentration is generally higher than those inferred from the 
CFHT observations, and more cons istent with those of the com- 
bine d weak and stro ng analyses of lUmetsu & Broadhursj ( 1200 8h 
and ISroadhurst et al . ( 2005b) , though not so high as the recent 
weak-lensing only ^edezinski et al, [( 120071) result, which all em- 
ploy the Subaru data. Overall, our results fit well into the body of 
weak lensing and combined strong & weak lensing work that ex- 
ists on the cluster. Compared to the strong lensing-only results, our 
spherical mass is low compared to those found using ACS data, as 
with all of the weak lensing results. 



Compared to the X-ray result, our mass and concentration are 
both high, though marginally consistent at la. This is unsurprising 
as i Ande rsson & Madejski (2004) find evidence for an ongoing or 
recent merger near t he cluster centre, a fin ding supported by strong 
lensing and flexion dLeonard et al.l ( l2007h ) measurements as well, 
and argue that X-ray analyses will underestimate the total mass in 
such cases. 



4.7.2 AbelI1835 

The catalogue employed for Abell 1835 is improved from that 
published in CS02 and C03, as the addition of two colour bands 
from CFHT observations allowed for a much more reliable dis- 
crimination between cluster, foreground, and background galax- 
ies. Using the catalogue constructed using a magnitude cut CS02 
reported NFW parameter values Af200 = 1.44 x 10^^ M© and 
C = 2.96. Using the new catalogue, and a lower mean source 
redshift < Zs >= 0.8, compared to the value of < Zs >= 1.0 
used in the CS02 and C03 analyses, we find for a spherical NFW 
marginalized mean parameters Af200 ~ 1.01±0.30x 10^^ Mq and 
C = 6.8lg I and approximate maximum-likelihood parameters 
Af200 = 'i-07to:ll X 10^^ Mq and C = 3.6l2;o- Our recovered 
concentration is higher, as expected given the significant improve- 
ment in the removal of faint galaxy contamination with the addition 
of two colour bands. The recovered mass, however, is lower, though 
still comfortably consistent with our new results at Icr. 

Looking to other weak lensing studies, our spherical mass 
is lower than that of Barde au et al. ( 2007), and our concentration 
higher. Howeve r , our mass and concentratio n are in agreemen t 
with the iDahld ( |2006|) result at Icr. Further, ISmith et al.l ( |2005|) 
find a total projected mass within 500 kpc of the cluster centre of 
M = (2.9 ± 0.6) X 10" /i-^Mq, for an Einstein-De Sitter cos- 
mology (they argue converting to a concordance cosmology would 
induce a change of no more than 10%, less than the error bars). For 
the same angular radius, our model in the concordance cosmology 
gives a projected mass of A/ = 3.0 X lO^"* /i-^Mq under a Flat 
prior, and for more direct comparison, Al — 2.7 x lO^** H^^Mq 
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under a Spherical prior, both in good agreement with the Smith re- 
sult. 

Our spherical mass and concentration estimates also agree 
very well with the parameter estimates derived from Chan- 
dra observations under a ssump tio ns of sphe r ical s ymmetry 
dSchmidt. Allen. & FabianI ( 1200 lb ). Izhang et al.l jlOOSh studied 
Abell 1835 with XMM and found a mass of M500 = (5.90 ± 
1.72) X 10" /i-^Mq (where Afsoo is the mass contained within 
a sphere with mean density 500 times the critical density); our 
marginalized mean model under a Flat prior gives M500 ~ 5.8 x 
10^* H^^Mq, also in good agreement. 



4.1.3 Abell 2204 

Though the colour selection in the catalogue used in this work is the 
same as that in CS02 and C03, as with the other two catalogues the 
PSF correction is improved and a lower mean background source 
redshift is employed to better reflect the expected mean redshift of 
the population after the colour selection. Our marginalized mean 
parameters (Afaoo = 0.69 ± 0.21 x 10^^ Mq- C = G.Gtll) and 
maximum-likelihood parameters (M200 = O-^'^to.ll ^ 10^^ M©; 
C = 4.5^2 4) under a Spherical prior agree with those of CS02 
at la. Our mass is lower than that reported by Dahle (20061), but 
with in their (large) la e rrors . 

IZhang et all ( I2OO8I) studied Abell 2204 with XMM and found 
an X-ray mass M500 = (4.12 ± 1.12) x lO" Ii^'^Mq- our 
marginalized model under a Flat prior gives iV/500 = 3.7 x 10^^ 
h~^MQ, fully consistent with the X-ray result. 

Thus, though there continues to be a significant scatter in the 
parameter values returned for all three clusters, our parameter es- 
timates under a spherical prior or in projection (where 3D struc- 
ture does not matter) are all comfortably situated within the errors 
of previous work utilizing weak lensing, strong lensing, and X-ray 
observations. The improvements in background galaxy discrimina- 
tion are apparent in the increased concentrations of Abell 1689 and 
Abell 1835. We are therefore confident our catalogues and MCMC 
sampler are returning reasonably accurate estimates of the mass 
profiles of these clusters, and understand how those profiles com- 
pare to those derived in previous work. Now we are equipped to 
examine the mass profiles of each of the clusters freed from the 
unphysical assumption of spherical symmetry. 

4.2 Triaxial parameter estimates 

The parameter estimates for halo mass and concentration using the 
triaxial NFW model should provide a better estimate of the true pa- 
rameters and errors of lensing galaxy clusters, as the model more 
closely reflects our true beliefs regarding the shape and structure of 
clusters. Generally, we find that triaxial models under the most gen- 
eral Flat prior return higher mass and concentration estimates than 
their spherical counterparts, but that otherwise the impact of triaxial 
fitting under other more stringent priors, whether drawn from the- 
ory or independent observational constraints, can shift parameters 
estimates in either direction from the typical spherical values. 

The question of which prior is most appropriate for a given 
question remains open, and crucially important for this undercon- 
strained problem. Generally, in the case where only weak lensing 
information is available, the most general Flat prior seems most 
useful, as it makes minimal assumptions about the halo shapes and 
parameter relations. Until the question of overconcentration is re- 
solved, it seems wise to avoid any prior that ties the concentration to 



the mass, as this may hide real discrepancies. Given that, the appli- 
cation even of a less controversial Mass prior is unwise, because the 
shape of the mass-concentration degeneracy automatically requires 
higher concentration values as mass values are pushed to their 
lower limits. This artificially exacerbates the mass-concentration 
relation discrepancy. A mixed AdC prior may be of great use in the 
future if the relation is better understood and constrained, as it can 
provide simple physical constraints on this underconstrained prob- 
lem without requiring independent observational data. However, it 
is only valid once it represents our true prior expectation for galaxy 
cluster structure, a status which current observational data does not 
yet confer. 

Similarly, employing a spherical prior imposes a prior con- 
straint that does not represent our true best predictions for cluster 
structures. The M200 — C contours show that generally the contours 
under the Spherical prior are enclosed by those of the Flat triax- 
ial models, though they occupy different subregions of the triaxial 
parameter space depending on the structure of the lensing cluster. 
Thus, they not only give errors that are too small and unrepresenta- 
tive of the true uncertainties in the problem, but unpredictably bias 
the mass and concentration estimates either high or low. Though 
these biases were shown to cancel with large-scale averaging in 
CK08, they are a large problem when studying individual clus- 
ters. Similarly, the overly-optimistic error estimates make proper 
comparisons of parameter values between different works inaccu- 
rate, as fully consistent results may appear excluded by inaccurately 
tight error bars; this is even more important because fitting spher- 
ical models to non-spherical data makes the details of the fitting 
method and its order of operations in measuring, averaging, and 
weighting shear far more important and likely to affect the final pa- 
rameter estimate. The inaccurate description of errors is also impor- 
tant when attempting to characterize the scatter in mass-observable 
relations, crucial knowledge for constraining the cluster mass func- 
tion and with it cosmological questions such as the normalization 
of the matter power spectrum and the nature of dark energy. 

As lOguri & Blandford ilOOj) convincingly demonstrated 
over a large population of triaxial halos, and as we argued in CK08, 
the appearance of circular symmetry in a projected cluster DOES 
NOT justify the assumption of sphericity, most especially in strong 
lensing clusters. Oguri & Blandford ( 2008) found that the strongest 
lensing clusters in the universe are expected to be more triaxial, but 
look more circular in projection, than their weaker counterparts, as 
line-of-sight projection boosts the lensing strength; CK08 showed 
that such halos leave no signature in their lensing signal to differ- 
entiate them from truly spherical halos, making it impossible to 
determine beforehand the degree of triaxiality of a given system. 
Thus, for all systems, and especially powerful strong lens systems 
such as Abell 1689, assuming spherical symmetry is never justified 
theoretically. 

When external information is available, such as strong lens- 
ing. X-ray, or dynamical data, this MCMC method allows for the 
straightforward combination of those external constraints with the 
weak lensing data. However, as emphasized in Sec. 13. II care must 
be taken that those external constraints do not reintroduce the very 
assumptions the triaxial model intends to avoid. In the case of 
strong lensing, this is most easily done through the 6 orientation 
prior, which comes from a robust theoretical prediction favouring 
line-of-sight oriented halos to have the highest cluster lensing sig- 
nals in the universe. The Einstein radius prior is in some ways 
preferable in that it requires no prediction from CDM theory, but 
is problematic in that, applied to a dark matter only profile such 
as the NFW, it neglects the impact of baryons in the very centre 
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Table 4. Mass and Concentration estimates from previous work. All lensing- and X-ray- derived values assumed spherical symmetry, and all X-ray analyses 
also assumed hydrostatic equilibrium. Values in parentheses were fixed during fitting. 



M200 [10^^ Mq] C Method Instrument Author 



A 1 f^9^C) 










1.47 


11 10 


WL 


ESO/MPG WFI 


CSOl method with current c3.t3.1ogu.e 


1.76 ± 0.20 


10.7+^-5+ 


WL 


Subaru 


Umetsu & Broadhurst (2008) 


- 


99 1+2.9 


WL 


Subaru 


Medezinski et al. (2007) 


-1 c:c-|-0.31 
1.00_Q 27 


4.28 ±0.82 


WL 


CFHT 


Bardeau et al. (2007) 


1.48 ±0.22 


7.6 ± 1.6 


WL 


CFHT 


Limousin et al. (2007) 


1 07+0-4,6 
^■'^'-0.36 


q C+0.5 

■^■"^-0.3 


WL 


CFHT 


Bardeau et al. (2005) 


1.13 


7.9 


WL 


ESO/MPG WFI 


C03 


1.07 


5.7 


WL (Infrared) 


NTT SOFI 


Kine et al. f2002b) 


85 


4 8 


WL 


FSn/MPG WFT 


Kino Plnwp Jt^ SrhnpiHpr Cyonj^ 


1.86 ± 0.16 


10.1+°| 

—0.7 


WL + SL 


Subaru + ACS 


Umetsu & Broadhurst (2008) 


2.25 ± 0.20 


7.6 ± 0.5 


WL + SL 


Subaru + ACS 


Halkola. Seitz. & Pannella (2006) 


1.72 ± 0.19 


10 9+1-^ 


WL + SL 


Subaru + ACS 


Broadhurst et al. (20(^ 




6.0 ± 0.6 (3tT) 


SL 


ACS 


Limousin et al. (2007) 


3.05 ± 0.30 


6.0 ±0.5 


SL 


ACS 


Halkola. Seitz, & Pannella (2006) 


3.19 


"■0-1.6 


SL 


ACS 


Broadhurst et al. (2005a) 


0.89«:- 


7 7+1.7 
'■'-2.6 


X-ray 


XMM 


Andersson & Madeiski (2004) 


A1835 












2.58 ± 0.48 


WL 


CFHT 


Bardeau et al. (2007) 


1.33 ±0.71 


(4.63) 


WL 


NOD/UHT 


Dahle (2006) 


1.44 


2.96 


WL 


ESO/MPG WFI 


CS02 


0.83 




X-ray 


XMM 


Maierowicz, Neumann & Reiprich (2002) 




■^■^'-'-0.63 


X-ray 


Chandra 


Schmidt, Allen, & Fabian (2001) 


A2204 










1.33 ± 0.82 


(5.03) 


WL 


NOD/UHT 


Dahle (2006) 


0.86 


6.3 


WL 


ESO/MPG WFI 


CS02 



of the cluster that are expected to significantly boost the strong 
lensing signal. This leads to a strong upward shift in the concen- 
tration value, which is unlikely to accurately represent the structure 
of the dark matter component of the cluster. X-ray constraints may 
be applied when a halo appears to have a low degree of triaxiality, 
making the assumption of spherical symmetry of the (always more 
spherical) 3D gravitational potential less problematic, and when the 
cluster appears relaxed and undisturbed. 

More detailed constraints may be applied by their direct inclu- 
sion in the likelihood function constraining the triaxial model. For 
example, strong lensing may be used to directly constrain the mass 
and concentration if the strongly lensed features, or the 2D conver- 
gence map constrained by them (constructed with no assumptions 
about the 3D structure of the halo), are used to directly constrain 
the 3D triaxial model. Doing so requires careful calibration of the 
weighting between the weak and strong lensing data, an optimiza- 
tion that has not yet been fully ex a mined but has been b egun by var- 
ious g roups ('Bradac et al] i l2005h ; ISradac et all ( l2006h ; iDiego et all 
( I2OO7I) ), and should be a focus of the next generation of galaxy 
cluster studies. 



4.3 Remaining systematics 

Though the parameters and the errors presented in this work now 
include the crucial contribution of halo triaxiality, there remain sev- 
eral outstanding systematics. 

One of these is uncertainty in the mean redshift of the back- 
ground galaxy population, changes in which would lead to changes 
in the mass estimates of several percent. Another is residual con- 



tamination by cluster galaxies; though the colour cuts made in all 
three catalogues are very stringent, there is still the possibility of 
contamination by a few cluster members. As with the uncertain- 
ties in the mean redshift of the background population, the effects 
should be negligible. 

More importantly, stringent colour cuts to remove contamina- 
tion often result in a significant decrease in the number density of 
background sources (for example, the catalogue from Abell 1689 
used in this work had a source density of only 7/arcmin'^). Low 
number density affects maximum-likelihood values and marginal- 
ized parameters significantly because of the large and non-Gaussian 
mass-sheet degeneracy between the NFW mass and concentration: 
mass and concentration constraints from low S/N data are poor be- 
cause the tails of the highly asymmetric M200 — C posterior prob- 
ability distribution elongate dramatically. Figure [TO] shows this ef- 
fect, taking 500 randomly chosen halves of the Abell 2204 cata- 
logue and treating them each as individual realizations to which a 
spherical NFW is fit. Red stars plot the distribution of recovered 
maximum-likelihood parameters, compared to the 2a confidence 
contour for the spherical model fit to the complete catalogue. The 
mean marginalized parameter values are even worse affected, with 
recovered values of concentration up to six times greater than the 
mean from the complete catalogue, because they are very sensi- 
tive to the long tails of the posterior when the signal is very weak, 
illustrated for five realizations in the Figure. The mean marginal- 
ized parameters would be expected to be more robust than their 
maximum-likelihood counterparts were the errors Gaussian; their 
sensitivity to the random removal of half the lensing signal reflects 
the very non-Gaussian errors induced by the mass-sheet degener- 
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200 ^ solar^ 

Figure 10. The red stars plot the maximum-hkehhood M200 and C pa- 
rameter values obtained fitting a spherical NFW to 500 randomly chosen 
halves of the Abell 2204 catalogue. The green circle gives the maximum- 
likelihood parameter values, and the shaded black region gives the 2a con- 
fidence region, obtained fitting the same spherical NFW to the whole cata- 
logue. The black contours give the 2a confidence contours for 5 of the 500 
half-catalogue realizations; the mass sheet degeneracy becomes very strong 
when the signal is weak, and causes the posterior probability distribution to 
develop very long tails to high concentrations and masses. 



acy. Because marginalized parameters are more meaningful in un- 
derconstrained problems such as that of fitting 3D mass models 
to lensing data, this is an important issue. More efficient ways of 
removing foreground contaminati on are therefore very valuable, 
su ch as the Bayesian met hod of [Limousin et al. I ( l2007h . or that 
of iMedezinski et al.l bOOTh that increases the efficiency of back- 
ground selection by studying the radial tangential shear profiles of 
various colour populations, or improved methods of colour discrim- 
ination based on detailed photometric studies of the background 
galaxy population. Future lensing work, including studies using the 
MCMC method employed here, can only be improved by increas- 
ing the background number density reliably available for inclusion 
in the lensing analysis. 

Our triaxial model includes no substructure, while we know 
galaxy cluster scale halos often contain subclumps, and certainly 
contain mass concentrations at the positions of the cluster galax- 
ies and their accompanying dark matter halos. The weak lensing 
data does not have the resolution o r signal-to-noise to constrain 
these; for example, iMarshalll ( l2006h found the Bayesian evidence 
supported fitting usually at most two components to cluster-scale 
substructures. In the future, in combination with strong lensing 
or flexion, where detail ed models including substructure are fre - 
quently constructed ( e.g. lLimousin et alJi2007h : |jullo et alj l l2007h : 
[Leonard et all j2007h ), it should be pos sible to constrain more of 
the s ubstructure of clusters. However, ' Clowe. De Lucia. & Kind 
( I2OO4) showed that the effects of substructure on parameter esti- 
mate are small compared to those of triaxiality, and so we expect 
the inclusion of multiple components to be a minor perturbation to 
the results we present in this paper. 

Finally, our quoted errors do not include a contribution from 



uncorrelated line-of-sight structure, shown by iHoekstrj ( l2003h to 
induce an unbiased scatter in parameter estimates that can in- 
crease the errors on concentration and mass by up to a factor of 
2. Dodelson ( 2004) demonstrated that this can be somewhat mit- 
igated by including a noise term for large scale structure in the 
lensing analysis; building on the greater coherence scale of large 
scale structure noise compared with the noise associated with the 
intrinsic ellipticity dispersion of galaxies, Dodelson notes that the 
errors on cluster mass can be reduced by ~ 50% for wide-field 
data, though this has yet to be implemented in any weak lensing 
analysis. 



4.4 Bayesian Evidence 

We choose not to use the Bayesian evidence to discriminate be- 
tween our fits under different priors, because we are not attempting 
to choose between them as models of the universe. The limitations 
of the evidence are clear when the Flat vs Spherical cases are con- 
sidered: in many cases the evidence would likely favour a Spherical 
over a Flat prior, due to a significant decrease in available parameter 
space coupled with a relatively small decrease in likelihood values 
under the Spherical prior, as seen in the very similar reduced 
values between all priors. However, unlike in cases of fitting multi- 
ple halos to account for substructure, or in another context, adding 
additional parameters to cosmological models, we know a priori 
from physical observations that a triaxial model is a better model 
than a spherical model - we see clearly in non-parametric lensing 
mass maps that galaxy clusters have significant levels of elliptic- 
ity and are not spherical! Thus, the prior on the spherical model is 
close to zero, but difficult to quantify. 

The various other theoretical priors studied in this paper are 
all derived from a single CDM model, and so treating them as sep- 
arate models to choose between is inappropriate. We study them 
here independently to understand how various aspects of the CDM 
model interact with the lensing data; in application to a large sur- 
vey, all priors that represent the true prior expectation for the prob- 
lem without interfering with the scientific questions asked by the 
study should be employed simultaneously. Similarly, external con- 
straints from strong lensing or X-ray do not represent different 
models to be selected between, only appropriate or inappropriate 
measures to constrain the triaxial models. 

The evidence is a powerful tool for model discrimination, but 
is not the appropriate method for choosing priors on a single phys- 
ical model, well-founded in theory and observations. 



4.5 Summary 

We fit a fully triaxial NFW model to weak lensing data from mas- 
sive galaxy clusters Abell 1689, Abell 1835, and Abell 2204, de- 
riving parameter estimates and errors under a range of theoretical 
and observational priors. Under a relatively weak prior drawn from 
the strong lensing orientation bias for Abell 1689, that same strong 
lensing prior combined with an X-ray derived mass and concentra- 
tion constraint for Abell 1 835, and under the very general Flat prior 
for Abell 2204, our best mean parameter estimates in a concordance 
cosmology with /i = 0.7 are 



• Abelll689: M200 

• Abell 1835: M200 

• Abell 2204: M200 



= (1.18±.23)xl0^^ Mq;C = 12.2±6.7 
= (0.98±.25) X lO^^^ (J ^ 3.7±i.O: 
= (0.72±.27) X lO^'^ M©; C = 7.1±6.2: 
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All results are consistent with previous work, but importantly with 
larger error contours that better reflect the true uncertainty in mass 
profile estimates. 

Triaxiality does not easily return Abell 1689 to the fold of 
average CDM clusters; the predicted CDM mass-concentration re- 
lation is enclosed in its 2a M200 — C confidence contour under 
most priors, but is excluded under some observational priors, mak- 
ing Abell 1689 at the moment weakly consistent with the predic- 
tions of CDM at 2(7. However, consistency does not mean truth: the 
primary result of this work is to demonstrate the necessity for im- 
proved methods of combining diverse observational constraints to 
counter the very large uncertainties that accompany fits of triaxial 
NFW profiles to galaxy cluster lenses. That way lies the future of 
galaxy cluster studies, whether of individual clusters or large popu- 
lations and mass-observable relations: without such optimized and 
robust combination methods, most results will continue to be con- 
sistent with one another and with theory simply through the size of 
their errors, rather than through physical concordance. 
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APPENDIX A: LENSING THROUGH TRIAXIAL HALOS 

Following OLS, the triaxial halo is projected onto the plane of the 
sky to find its projected elliptical isodensity contours as a function 
of the halo's axis ratios and orientation angles (6^, (f}) with respect 
to the the observer's line-of-sightQ The elliptical radius is given by 
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where (X, Y) are physical coordinates on the sky with respect to 
the centre of the halo. 
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The axis ratio q of the elliptical contours is then given by 

1=^ (A8) 

qx 

and their orientation angle \& on the sky by 

* = i tan-' -^-^ {qx > qv). (A9) 

Here we diverge slightly from OLS's treatment as we are in- 
terested not in deflection angles but in the lensing shear and con- 
vergence, both combinations of second derivatives of the lensing 
potential <E> (commas indicate differentiation): 
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(All) 
(A 12) 



These derivatives are calculated as functions of integrals of the 
spherical convergence (see e.g. Bartelmann ( 1996) for a full 
treatment of weak lensing by a sph erical NFW pr ofile) following 
the method of lSchramnj ll990l) and lKeetoiJ j200lh . normalised by 



a factor of l/\/7 from Equation I A4I (see OLS for the derivation of 
this normalisation) 
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Note that our radial variable C, appears different from Keeton's 
5 because it is defined in terms of two axis ratios qx and qy rather 
than one q: C, = ^/qx- This reflects a dependence on the 3D struc- 
ture of the cluster; for example, extended structure along the line of 
sight decreases qx and thus increases the convergence and shear at 
a given {X, Y). 



^ Although we set c = 1, we keep c as a variable in our notation for 
consistency with OLS. 
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APPENDIX B: PRIORS FROM SIMULATIONS: SHAW PRIOR 



We define the Shaw prior, plotted in Figure[T] by fitting polynomials to the data points of Figure 14 in lShaw et alj j2006h : 

f if 6 < 0.5 

^ I [l-6329fe^ - 7.97756* + 9.3414fo^ - 6.6558fe^ + 2.29646 - .3088] x 10^ if 0.5 ^ 6 ^ 1.0 
and 



^ J 5.76647 (f)®- 2.459265 (f)^+ 42.3154 (1)"- 37.2765 (f)^ if 0.65 s: f ^ 1.0. 
17.4650 (f )^ - 4.00238 (f ) + .32462 x 10"* 



if f < 0.65 
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